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Abstract 

Quantum Brownian motion in the strong friction limit is studied based on the exact path integral 
formulation of dissipative systems. In this limit the time-nonlocal reduced dynamics can be cast into 
an effective equation of motion, the quantum Smoluchowski equation. For strongly condensed phase 
environments it plays a similar role as master equations in the weak coupling range. Applications 
for chemical, mesoscopic, and soft matter systems are discussed and reveal the substantial role of 
quantum fluctuations. 



PACS numbers: 03.65.Yz, 05.40.-a,73.23.-b,82.20.-w 



Quantum Brownian motion is much more involved than its classical analog 
since in general tractable equations of motion do not exist. Progress has been 
made in the weak friction regime leading to a variety of master equations. Here, 
we analyze the opposite domain of very strong friction and reveal that a gen- 
eralization of the classical Smoluchowski equation, the quantum Smoluchowski 
equation, can be derived from first principles. This opens the door to study 
the dynamics of strongly condensed phase systems at lower temperatures. In 
particular, quantum fluctuations turn out to play a substantial role as shown 
explicitly for three examples from chemical reactions, mesoscopic physics, and 
charge transfer in molecules. 



I. INTRODUCTION 



Brownian motion, that is the fate of a heavy particle immersed in a fluid of lighter 
particles, is the prototype of a dissipative system coupled to a thermal bath with infinitely 
many degrees of freedom. The work by Einstein in 1905 has developed a mathematical 
language to describe the random motion of the particle and has uncovered the fundamental 
relation between friction, diffusion and the temperature T of the bath. Half a century 
later, this seed had grown into the theory of irreversible thermodynamics '2\, which governs 
the relaxation and fluctuations of classical systems near equilibrium. At that time, a new 
challenge had emerged, the quantum mechanical description of dissipative systems. 

In contrast to classical Brownian motion, where right from the beginning the work by 
Einstein and Smoluchowski had allowed to consider both weak and strong friction, 

for a long time, the quantum mechanical theory could handle the limit of weak dissipation 
only. In this case the interaction between the "particle" and the "bath" can be treated 
perturbatively and one can derive a master equation for the reduced density matrix of the 
"particle" ^. This approach, based on the so-called Born and Markov approximations, has 
been very successful in quite a number of fields emerging in the fifties and sixties of the last 
century, such as nuclear magnetic resonance ^,0] and quantum optics j?]. 

Roughly, a dissipative quantum system can be characterized by three typical energy 
scales, an excitation energy hujQ, where Uq is a characteristic frequency of the system, a 
coupling energy hj to the bath, where 7 is a typical damping constant, and the thermal 
energy fc^T. The weak coupling master equation is limited to the region h'y <^ huo^, ksT. 
This is the case whenever the typical linewidth caused by environmental interactions is small 



compared to the line separation and the thermal "Matsubara" frequency 2TikBT /h. On the 
other hand, the approach will fail for strong damping and/or low temperatures. 

The work by Feynman and Vernon P] has shown how to take advantage of the path 
integral representation of quantum mechanics to derive an expression for the reduced density 
matrix of the "particle" as a sum over forward and backward paths, which is valid for 
arbitrary damping strength and temperature. These two sets of paths arise from the two 
time evolution operators in the formal expression for the time-dependent density matrix 
W{t) = exp[—{i/h)Ht]W{0) exp[+{i/h)Ht], where H is the Hamiltonian of the system. 
In the path integral representation of the reduced density matrix the influence of the bath 
emerges as time-nonlocal terms in the action governing the path probability. These terms 
include also a coupling between the forward and backward paths and make it difficult to 



evaluate the path integral explicitly, although progress can be made in some cases 

Everything becomes simpler in the limits of weak and strong damping. While the former 
case leads to the master equation approach mentioned above, the latter case is usually 
referred to as the Smoluchowski limit. In the classical case this limit is well understood and 
simplifies matter considerably, since the momentum P of a heavily damped particle with 
position Q is a slow variable that can be eliminated adiabatically. It is intuitively clear, 
that a quantum mechanical theory cannot entirely dispose of the variable P conjugate to Q. 
In fact, for a classical system with characteristic frequency ujq the roots X± = ±iu!o of the 
characteristic equation in the undamped case are in the presence of damping turned into 
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For large friction, j ':$> ujq, this can be approximated by 

A± = 




While in the strong friction limit the slow classical dynamics is governed by the small root 
c<Jq/7, quantum mechanics accompanies reduced fiuctuation of Q with enlarged fiuctua- 
tions of the conjugate variable P. Both roots X± are essential to determine the fiuctua- 
tion spectrum. That is why quantum effects are not restricted to the temperature range 
ksT < huj^/'-f, rather, in quantities sensitive to fluctuations, quantum effect can manifest 
themselves in the strong friction limit even up to higher temperatures than in the weak 
damping limit. This makes the quantum Smoluchowski limit addressed in this article non- 
trivial. 



The paper is organized as follows. In the next section|n]we shortly review the formulation 
of dissipation within the path integral formalism and sketch the derivation of the Quantum 



Smoluchowski Equation (QSE) [11|. The second part (Sec. Illl|) is devoted to three specific 
applications. 

II. THE QUANTUM SMOLUCHOWKSI EQUATION 
A. Reduced dynamics 

As already pointed out in the Introduction the inclusion of quantum dissipation in a 
non-perturbative way has been established only since the early 80s. The standard approach 



12| | starts from a system+reservoir model 

H = Hs + Hji + Hj (1) 

with a system part Hs, an environmental part Hji, and a system-bath interaction Hj. The 
reservoir (heat bath) is mimicked by a quasi-continuum of harmonic oscillators bilinearly 
coupled to the system. Dissipation appears when one considers the reduced dynamics by 
properly eliminating the bath degrees of freedom. In fact, one this way regains in the 
classical limit a generalized Langevin equation. It is thus important to realize that the 
only restriction associated with the oscillator bath is the Gaussian stochastic nature of the 
environment, which in turn means that the relevant impact of the bath onto the system 
dynamics is completely determined by the average and the autocorrelation function of the 
bath force. 

The quantum dynamics of the reduced system follows from 

p{t) =TTR{exp{-iHt/h) W{0) exp{iHt/h)} (2) 

where W{0) describes the initial state of the total system. In the ordinary Feynman- Vernon 
theory this state is assumed to be a factorizing state, 1^(0) = ps{0) exp{—f3Hfi)/Zji {Zji 
is the bath partition function and (3 = 1/k-QT), so that each one, system and equilibrated 
bath, lives in splendid isolation at t = 0. While this assumption may be justified in the weak 
damping/high temperature limit, it certainly fails for moderate to strong friction and/or 
lower temperatures. It can be shown explicitly that in the classical limit even the Langevin 
equation is not regained, but differs by initial boundary terms that may persist up to long 
times. In the strong damping realm considered here, the initial state must be correlated 



[10|, i.e. W{0) = J2j Oj exp{—j3H)0j/Z with preparation operators Oj and Oj acting onto 
the system degree of freedom only and the total partition function Z. For transparency, in 
the sequel we assume these operators to be diagonal in position so that the initial density 
for a one-dimensional system reads in position representation 

PiQi, Qi, t = 0) = ppiQi, q[) X{qi, q[) (3) 

with the reduced thermal equilibrium density matrix Pp^q.q') = Ti R{q\exp{— (3 H)\q') and 
the coordinate representation X{q,q') = J2j{(l\0j\q) {q'\Oj\q') of the preparation operators. 
For a specific initial state, one does not need to know the specific form of these operators 
(which may be very complicated), but rather chooses the preparation function properly. 

The reduced quantum dynamics starting with an initial state is now obtained 
within the position representation by employing the path integral formalism. Since the bath 
contains harmonic degrees of freedom only, it can be integrated out exactly and one arrives 
at 

pilf^Qf^t) = J dqidq- Jiqf,q'f,t,qi,q'i)X{qi,q'i} (4) 

where the propagating function J(-) is a threefold path integral — two in real time, one in 
imaginary time — over the system degree of freedom only 

J{qf, q'f, t, g., Q'i) = y mW]m e^^t'^''^' '^1/^ (5) 

with Z = TT{exp{~(3H)}/Zji. The two real time paths q{s) and q'{s) connect in time t the 
initial points qi and q'^ with fixed end points g/ and q'j, while the imaginary time path q{a) 
runs from q^ to g- in the interval hf3. The contribution of each path is weighted with an 
effective action S[g, q', q] = Ss[q] — Ss[q'] + iS[q] + i(t>[q, q,' , q] which consists of the actions 
of the bare system in real and imaginary time, respectively, and an additional interaction 
contribution (influence functional) non-local in time. The latter one can be written as 

(f)[q] = fdz f dz' q{z) K{z - z') q{z') + '-p j dzq{zf (6) 

where the ordered time integration is understood along the contour: z = s from t — s> 0, 

—iTij3 + s from ^ t with 

q'{s) for z = s < s <t 

< g(r) for z = -IT < t < hf3 ■ (7) 

q{s) for z = —ihp + s < s <t 



z = —IT from — >• Tifi, z 



q{z) 



The effective impact of the bath is completely controlled by the damping kernel 

^ ' Jo 71 ^ ' sinh{ujh(3 /2) ^ ' 

where I{uj) denotes the spectral density of the environment. As expected hK{z) coincides 
with the autocorrelation function of the bath force. In particular, for real times the kernel 
K{s) = K'{s) + iK"{s) is related to the macroscopic damping kernel entering the classical 
generalized Langevin equation 

2 duj liuj) , , , , 
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via K"{s) = {M/2)d'y{s)/ds and K'{s) M'y{s)/h(3 in the classical limit (M is the mass 
of the Brownian particle). The term with /i = \imfii3_Qhj3K{0) in (jH)) gives a potential 
renormalization due to shifts of the minima of the bath oscillators by the coupling to the 
system. 



B. Quantum Smoluchowski range and thermal equilibrium 

The classical Smoluchowski limit is related to a separation of time scales between fast 
equilibration of momentum and slow equilibration of position. In this way, the Fokker-Planck 
equation for the phase space distribution can be adiabatically reduced to a Smoluchowski 
equation for the marginal distribution in position space [13]. For quantum dissipative sys- 
tems the expectation is that friction makes the system to behave more classically so that for 
strong friction the complicated path integral expression simplifies considerably. This is 
indeed the case as we have shown recently ^3] • In the sequel we will briefly collect the main 
findings and, to keep things as transparent as possible, consider a one- dimensional system 
only moving in a sufficiently smooth potential field V{q). 

A typical damping strength in the long time limit is defined as 

where 7(ci;) is the Laplace transform of the classical damping kernel 7(s). For ohmic friction, 
I{uj) = M'yuj, for instance, one finds 7 = 7. The same is true for the more realistic Drude 
damping /(cu) = M^uuj'^/^uj'^+u!'^) with cut-off frequency Uc- Now, given a typical frequency 
uo of the bare system (e.g. its ground state frequency) by strong damping we then mean 

(cf. fig. CD 
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FIG. 1: Quantum Smoluchowski range according to (jllj) (shaded). The classical range 'jhp < 1 is 
simple shaded, the quantum range 'yh(3 > 1 double shaded. 

Hence, we extend the time scale separation known from the classical Smoluchowski range to 
incorporate the time scale for quantum fluctuations hjS. As we will see, this does not mean, 
however, that quantum effects are negligible at all, since they are related to Matsubara 
frequencies z/„ = 2Tm/h(3, = 1, 2, 3, . . . which can be arbitrarily large. 

The idea is now, to evaluate for strong friction 7/^70 ^ 1 the path integral expression 
(HJ on a coarse-grained time scale s ^ Ti[3, l/c<Jc, I/7 and r ^ l/cUc, I/7. The consequences 
are the following: (i) Non-diagonal elements of the reduced density matrix are strongly 
suppressed during the time evolution, (ii) the real time part of the kernel K{s) becomes 
local on the coarse grained time scale, and (iii) initial correlations described by (jHl) survive 
for times of order 7/ix'o verifying that factorizing initial states cannot be used. 

In a first step, this program is applied to calculate the thermal (unnormalized) equilibrium 
density matrix 



with E[g] = E[0, 0, q\. Due to (i) the strategy is to invoke a semiclassical type of approxima- 
tion by assuming self-consistently that |^(t) — ^(0) | remains small on the time interval hl3. In 
function space the path integral is then dominated by the contributions of the minimal action 
paths and Gaussian fluctuations around them. The former ones are determined perturba- 
tively up to corrections of order I/7 and its corresponding minimal action — conveniently 
expressed in terms of the sum f = {q + q')/2 and difference coordinate x = q — q' — follows 
as 




(12) 
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with 
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(14) 



Note that in A the quadratic dependence on i/„ originating from the inertia part in the 
equation of motion and in Vt the cut-off frequency ujc appearing in ^{uj) are essential for 
the convergence of the corresponding sums. This reveals a fundamental difference to the 
classical Smoluchowski limit and indicates the importance of quantum fluctuations even for 
strong friction. From the above, one observes that A measures the typical strength of quan- 
tum fluctuations in position and VL is related (via the identification %/ Ax — > Ap) with the 
equilibrium variance in momentum Vt = (p^). Before we discuss further details, let us arrive 
at the complete expression of the thermal equilibrium density matrix. The remaining Gaus- 
sian fluctuations around the minimal action paths are treated by switching to Matsubara 
frequency space. Correspondingly, a straightforward calculation provides together with 



where Z denotes a proper normalization factor, e.g. Z = J dqpp{0,q). This is the first 
important result of the strong friction analysis: The equilibrium distribution of a strongly 
damped quantum system in arbitrary (sufficiently smooth) potentials. In particular, for a 
harmonic oscillator the known result is recovered To be specific we consider Drude 

damping j{z) = ^ujJ [oJc + z) with a high frequency cut-off ^ 7- Then, the functions A 
and can be expressed in terms of \E' functions 



where C = 0.577. . . is Euler's constant. In the high temperature limit ^yhf] <^ 1, one finds 
A ^ h'^P/12M and ^ M/P so that the Wigner transform of reduces to the classical 
phase space distribution. The dependence on friction appears for lower temperatures as 
a genuine quantum effect. Of particular interest is the limit 7^/5 ^ 1 where we have 
A ^ {h/Mj'7i)\og{'yh/3/27!-) and Q ^ {Mh'y / Tr)\og{uJc/ ■ Quantum fluctuations in position 
are suppressed by friction, not algebraically though, but much weaker. They also show 
a nonlinear dependence on h meaning that for 'jh(3 ^ 1 we work nevertheless in a deep 




-l3V(f)-nx'^/2n'^ A/3[/3V"(f)2/2-3V"(f)/2] 



(15) 




(16) 



quantum domain. The same can be seen from the momentum variance which grows with 
friction to guarantee Heisenberg's uncertainty relation. Accordingly, as assumed above, off- 
diagonal elements of the distribution pp{q, q') are strongly suppressed with \x\ being of order 
1 / ^J^y\og{LjJ^ or smaller. Concluding this section, we mention that higher order corrections 
in A are associated with higher than second order derivatives of the potential V{q). 



C. Quantum Smoluchowski equation 

The analysis of the previous section already indicated the strong suppression of off- 
diagonal elements of the density distribution in thermal equilibrium. While in principle 
we could now proceed to study the reduced dynamics in real time for the total distribu- 
tion p{qf,q'j-,t) y^], we concentrate in the sequel on the quantum analog of the classical 
Smoluchowski limit and thus restrict ourselves to the time evolution of its diagonal part, 
the position distribution P{qf,t) = p{qf,qf,t) In fact, one can show that off-diagonal 
elements relax to thermal equilibrium with respect to the instantaneous position of the 



Brownian particle on a time scale of order I/7 15|. For this purpose it is convenient to 



introduce sum and difference paths also in real time 

r{s) = [q{s) + q'is)]/2 , x{s) = g(s) - q'{s) (17) 

with corresponding fixed end-coordinates Vf = r{t) = qf,Xf = x{t) = and initial coordi- 
nates rj = r(0),Xj = x(0) distributed according to the initial state Q. Already the classical 
Smoluchowski limit reduces the set of acceptable initial states to those with momenta suffi- 
ciently bounded from above. In the quantum domain, the detailed analysis reveals that the 
Xj-dependence of the preparation function A(xj,rj) must be sufficiently smooth within the 



range \xi/yh/M'y\^^hp/{'y/ujQ) such that we may put A(xj,rj) ^ A(0,rj). 

Now, let us first look at the leading order contribution where quantum fluctuations in 
position are neglected, i.e. A = 0. This is expected to lead us back to the classical Smolu- 
chowski equation, for 'jhjS ^ 1 though, i.e. far from the classical limit 7^/5 <^ 1. As above, 
we assume excursions from diagonality, encoded in the x-paths, to be small, roughly at 
most of order \ j By expanding the effective action a;,r] up to second order in x we 
arrive at a solvable Gaussian path integral for the x-paths. We put x(s) = + (5x(s) with 
deviations 5x(0) = bxit) = and obtain to leading order 



with the action 



Sir] = f ds [M^r + V'{r)f . (19) 
Jo 



Thus, the imaginary time path contribution follows from the leading order expression of 
fll5|) for X = Xj = 0, Tj = f. One is left with the r-path integral which, written in form of a 
propagator, appears as 

P(g/, t) = J dq, G{qf, t, q^) P{q,, 0) (20) 

where in G{-) the paths r{s) run from r(0) = qi to r{t) = qf, each one weighted with the 
exponential on the right hand side of ()18|) . One this way ends up with the path integral 
representation of the classical Smoluchowski dynamics where P{q, t) obeys 

dPii,t) _ 1 d 

^^-M^a^^^'^^^''^ ^^^^ 
with the classical Smoluchowski operator Ld = V + k-oTdq. Indeed, from this derivation 
one finds x{sY to be at most of order h^P/M'yt. 

To include quantum fluctuations to lowest order in A, we could proceed to directly cal- 
culate corrections to the propagator as in 11], which is rather involved, however. Here, we 
follow a somewhat simpler route and first exploit that the flux vanishes in thermal equilib- 
rium. Accordingly, one writes dfP = {l/'yM)dq{l + 5f)Lq^P with a dynamical correction 
(5/, and an operator Lqm determined perturbatively up to order A by LqmP/j = O(A^) with 
the equilibrium density specified in fll5|) . The latter one already contains the dominant part 
of the quantum fluctuations in the QSE (see below). To verify that 5f can be disregarded, 
it suffices to study the path integral dynamics in the time range h(3, 1/7, l/c^c ^ t <^ 1 /oJq 
where it becomes effectively Markovian. Then, along the lines described above a semiclas- 
sical approximation applies and leads to a propagator including A dependent corrections. 
From this, one finds 5f oc 0(A/7^) and thus to be negligible. As the result, the equation 
for P(g, t) is of the form of the classical Smoluchowski equation ()21|) but with L^i replaced 
by -£qm, i-e., 

Kr. = UM + k^Tdqb{q) (22) 
where we introduced an effective force field and a diffusion coefficient, respectively, 

U^q) = V'{q) + ^V"{q) , b{q) = 1 + ApV"{q) . (23) 

The expressions and (^Hj) determine the dynamics of overdamped quantum systems 
at lower temperatures. The corresponding quantum Smoluchowski equation has still one 
deficiency though, which is directly related to the fact that it has been derived perturbatively. 



Namely, it conserves thermodynamic symmetries, as e.g. vanishing stationary currents in 
thermal equilibrium, only to order A meaning that in a direct numerical integration higher 
order terms may give rise to unphysical results. The most obvious way to see this, is to 
consider the stationary current J^t = limt^oo{<i) /d in a periodic, asymmetric potential V{q + 
a) = V{q),V{q) ^ V{—q) with no external bias (ratchet potential) It is determined 

from the stationary distribution Pst via LqmPst = —iMJ^t as 
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dq 

D{q) 
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dy e 



my) 



(24) 



where ipi^q) = Jq dyU^fi{y) / D{y) . Apparently, the diffusion coefficient in leads to ip{a) = 
O(A^) and thus to Jst 7^ 0. The strategy to cure this problem is to look for a "uniformization" 
of the operator Lqm by finding a proper diffusion coefficient D{q) which coincides with D 
in order A, leads to the equilibrium distribution (P3|) . and, of course, respects Jgt = in 
thermal equilibrium. The only functional form consistent with these constraints brings us 
to 

D{q) = l/[l-Af3V"{q)] . (25) 

We note that D can be seen as a Pade approximant to a series in |A/?V^"| < 1 starting with 
the terms included in D. The corresponding expression ()24|) for the stationary current is 
now given by the potential ipiq) = ip{q) — with 



(26) 



which obeys full periodicity in agreement with the correct symmetry Jgt = 0. Further, 
Pfi{q) oc exp[— /5'?/'(g)]/D(g) coincides with p/3(0, q) (fT3j) up to corrections of order A^ in the 
exponential which is consistent with our perturbation theory. 

Eventually, we arrive at the quantum generalization of the classical Smoluchowski equa- 
tion, the so-called quantum Smoluchowski Equation (QSE), 



dP{q,t) 



1 d 



UM + k^T^D{q) 
oq 



(27) 



dt -fM dq 

with Ucs specified in and D in (j^Hl)- The quantum analog to the classical Langevin 
equation in the strong damping limit follows (in the Ito sense) as 



M^q + U',M=atWDiq) (28) 

with Gaussian white noise (^(t)) = 0, {^{t)^{t')) = 2M-fk^T5{t - t'). Eqs. (EB and ^ 
are the main results of this study. They describe the reduced dynamics of an overdamped 



quantum system from high down to very low temperatures and show that the corresponding 
quantum stochastic process is equivalent to a classical process in an effective potential and 
with multiplicative noise. 



III. APPLICATIONS OF THE QUANTUM SMOLUCHOWSKI EQUATION 

The QSE plays in the strong friction limit a similar role as the quantum master equa- 
tions in the weak damping limit, since it allows for explicit results in strongly condensed 
phase systems. Three examples from chemistry, mesoscopic physics, and soft matter will be 
presented in this second part. 



A. Escape from a Metastable Well 

The problem of escape from a metastable well can be found in a variety of realizations 
including chemical reactions, diffusion in solids, or nuclear fission processes, to name but 
a few (see also the contribution by Hanggi and Ingold Its particular feature is a 

separation of time scales between local well motion and long time decay characterized by a 
decay rate. Classically, a first thorough analysis has been performed in a seminal work by 
Kramers with substantial extensions since then For quantum mechanical systems a 
large amount of research has been done in the 1980s in the context of macroscopic quantum 
tunneling and the field has gained renewed interest recently due to its importance for solid- 
state based implementations of qubits. Here, we look at the problem at strong friction and 
low temperatures. 

The archetypical situation is the following: In a metastable potential V{q) a barrier 
with height Vb ^ k-gT, Kuoq separates a well region (well frequency uoq) from a continuum. 
Initially particles stay in local thermal equilibrium inside the well. As time elapses, particles 
surmount the barrier and for intermediate times (plateau range) their position distribution 
becomes quasi-stationary P{q,t) —>■ Pst(Q') describing a constant flux across the barrier 

<^st = ~^^^^-^qm-Pst • (29) 

The decay rate follows from 

r = (30) 

■^well 

with the well population Z^^n. For strong friction the changeover from quasi-equilibrium in 
the well around q = and nonequilibrium on the other side of the barrier is restricted to 
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FIG. 2: Ratio of the exact rate with the classical (dashed) and the QSE (solid) rate vs. friction 
for \V"{qb)\ = V"{0) = MujI 

the vicinity of the barrier top located at qb > 0. Hence, the stationary distribution takes the 
form Pstiq) = Pfi{(l) 9st{(l) with a form factor obeying Qstiq) ^ 1 in a close range to the left 
of qh and gst{q) ^ in a close range to the right of g^. In fact, from (|^ a direct integration 
yields 



with if) specified in fl2()|l . We note in passing that this distribution is identical to that 
obtained from real-time calculations [ly]. Now, inserting Z^^w = I-oo'^1-Pstic[) into (jHU)) 
and by invoking a harmonic approximation around the well minimum and the barrier top 
according to the above discussion, we obtain the rate in the quantum Smoluchowski range 

as 



In this expression the second exponential accounts for quantum fluctuations, while the first 
factors coincide with the overdamped Kramers rate. Note that A dependent terms enter 
exponentially and thus substantially enhance the quantum rate compared to the classical one 
(=ee fig. a. Part.culariy in the high temperature domam the tate enhancement . damping 
independent and takes the form already derived by Wigner |l2j. In the low temperature 
domain, however, a complicated dependence on damping appears associated with a strong 
increase of the rate compared to its classical value. A quantum rate expression that becomes 
exact in the semiclassical limit of a high barrier has been derived in [2i| . For strong friction 
this formula reduces to ()H2|1 and already for moderate friction agrees well with the QSE- 




(31) 




(32) 



result, fig. |21 Remarkably, the rate enhancement is observable already at relatively high 
temperatures ^V"{0)/Mhl3 < 1 provided damping is sufficiently strong to guarantee •yhp 2> 
1. 



B. Phase Diffusion and Charging Effects in Josephson Junctions 

Since the discovery of the Josephson effect, devices based on Josephson junctions (J J) have 
revealed an extraordinary wealth of phenomena studied theoretically and experimentally 



as well 



21 



l22l |23| . Recent realizations of the underlying system, two superconducting 



domains separated by a tunnel barrier, comprise superconducting atomic contacts j2^ and 
solid-state based quantum bits 2^. In essence, two parameters determine the dynamics, 
namely, the coupling energy of the adjacent domains (Josephson energy) Ej and the charging 
energy of the contact Ec = 2e^/C of a junction with capacitance C. It turns out that 
the competition between these two scales is crucially influenced by the electromagnetic 
environment surrounding the junction, i.e. its impedance which in the simplest case is given 
by an ohmic resistor with resistance R. 

In the context considered here, a particular feature of a JJ is the fact jjiat its dynamics 
can be visualized as a diffusive motion of a fictitious particle (RSJ model) |2l|. This mapping 
is due to the famous Josephson relations 

/s = /csin(0) , 0(t) = |l^(t) (33) 

with Ic = {2e/h)Ej, the difference between left and right superconducting phases 0, and 
where the phase velocity is related to the voltage drop V{t) across the junction. Charge 
transfer captured by these relations corresponds to a phase coherent Cooper-pair current. 
With the translation rules to the mechanical analog 

the diffusive dynamics of the phase is determined by a generalized Langevin equation. The 
case of very small capacitance corresponds to very strong friction so that for a current biased 
JJ the classical time evolution reads 

M70 + dU{(j))/d(j) = ^{t) (35) 



with the Josephson potential U{(j)) = —Ej cos(0) — Ei(f) where the energy Ei is related to 
the bias current {h/2e)I, and current noise (^(t)) = 0, (^(t)^(tO) = i'^l/P)^^ - t')- 



This classical Smoluchowski dynamics has been studied in detail already in the 60s 
The QSE developed above, provides a generalization of this description to lower 



temperatures where charging effects enter the game 28[. In fact, it has been shown in 
the context of single charge tunneling that for sufficiently low temperatures these lead to 
Coulomb blockade associated with incoherent transfer of Cooper pairs through the junction 

We start by translating the constraints for the QSE dynamics to the case of JJs. Here, 
in addition to the relation we also have to take into account that typically the junction 
is subject to an external voltage (or a corresponding current). Accordingly, in order for 
the momentum M0 to relax within the RC = I/7— time to a Boltzmann-like distribution 
around (0), the external voltage V is restricted by eV <^ hj. By combining 'j'^/uJq ^ 1, 7^/3, 



where uq = \^2E^E^/h is the plasma frequency of the unbiased J J, with this latter condition 
and expressing them in junction parameters we expect the QSE to capture quantum phase 
diffusion in JJ if 



Ej27rV ' 27r2p' RI^ 
In the above, p = R/Rq with the resistance quantum Rq = h/Ae^. Since typically p ^ 1 the 

above condition allows for a broad range of values for Ec/E;}, and also large voltages 

V/RIc- We note in passing that the above relation also ensures that the actual non-ohmic 

impedance seen by the junction can effectively be treated as ohmic. 

Now, with the Josephson relations (jHHj) the effective potential and the corresponding 

diffusion constant entering the QSE (j2Z|) are found as 

t/eff(0) = -E*j cos(0) - Ei0 , = Ej[l~ I) (37) 

and 

D(0) = [l-e cos(0)]~^ , e = A/5Ej . (38) 
From the crucial A-function specified in and reading here 



A = 2p 



27i'p , / PE, 



(39) 



PE, \2ti^p^ 

one sees immediately that (3Ec/ttp = 'yhf3 controls the changeover from classical to quantum 
Smoluchowski dynamics. In particular, it turns out that barrier related quantum fluctuations 
incorporated in E^ dominate in the limit of jSEj <^ 1, while diffusion related ones, accounted 
for by the parameter 6, prevail in the opposite range ^ 1. 

The important observable is now the response of the JJ to an external bias cur- 
rent, namely, the average voltage (V). It is related to the steady state current via 
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FIG. 3: Range of the QSE for a J J with p <C 1, V/RIc < 1. The classical IZT range (shaded) and 
the domains of Coulomb blockade (CB) and macroscopic quantum tunneling (MQT) are indicated. 
The QSE is applicable above the thick line, see (j36|) . and the arrows illustrate various changeovers 
discussed in the text. 

{V) = {h/2n)\imt^oo{<P) = 27rJst- Jst reads as in with periodicity a = 2it, the dif- 
fusion constant from p8|l . and the potential (j22|l adapted according to (j37jl and p8|l . One 
then finds the current voltage-characteristics of a current biased junction to read 



(40) 



Here, the nominator Tq^ results from normalizing the steady state phase distribution to 1 
and can be written as 



27r JO Jo 
x[l-e sin(0' - 0/2)] e2^^«('^'<^') 



(41) 



with 



^(0, 0') = sin(0') sin(0/2) [Ei + E^j cos(0') cos(0/2)] 



(42) 



The expression ()4()|1 together with ()4H1 is the central result from which various known 
findings can be derived as limiting cases (see fig. Ej). (i) For (3Ec/ p ^ 1 the function Tq^ 
reduces to its classical form (i?j Ej, 6—^0) and the classical Ivanchenko-Zil'berman 
Theory (IZT) Q is recovered, (li) In the low temperature domain jSEc/p ^ 1, but for 
smaller couplings jSEj < 1, we have 6' ^ 1 so that diffusion related quantum fluctuations 
are negligible. As a consequence, quantum effects in the transport can be described by a 




FIG. 4: Left: Current-voltage characteristics for PEj = 2, p = 0.04. The quantum results [PE^. = 1 
(dashed), PEc = 20 (solid)] are shown with the classical one {PEc = 0, dotted). Right: Super- 
current vs. voltage for PEj = 0.25, p = 0.04. The classical result (dotted, from [2a|) and the 
CB expression [dashed, from 31]] are depicted together with the QSE result [solid, from ^]] for 
PEc = 0.15 (top), pEc = 20 (middle), PE^ = 1000 (bottom). 

renormalized coupling energy Ej. This important extension of IZT has first been derived in 
31[ based on a direct evaluation of the real-time path integral expression. The supercurrent 
across the junction coincides with results from Coulomb blockade (CB) theory (cf. fig. 
thus describing an incoherent transfer of charges, (iii) For PEc/ p ^ 1 and sufficiently larger 
couplings PEj > 1, coherent Cooper pair tunneling exists. Then, for a = I / Ic < 1 occasional 
phase slips occur and lead to the voltage 



{V) \/l - „-2/3£;j(l-a2)3/2/(3„2) 29Vr^ 
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which via 6 is affected by diffusion related quantum fluctuations. As can be observed in 
fig. El the result PHj) tends for ^ to classical thermal activation over the barriers 
of the washboard pot ential f/(0), where quantum corrections are of the known damping 
independent form [S^. At lower temperatures, i.e. for finite 9, they show a complicated 
dependence on p and capture the precursors of macroscopic quantum tunneling (MQT) 
found at very low temperatures 33|. Thus, the central result pn|) fills the gap between 
established results in different transport domains: On the one hand, for fixed PE^/p > 1 it 
leads with increasing PEj from Coulomb blockade to coherent Cooper pair tunneling [fig. 13 
arrow (a)]. On the other hand, for fixed PEj > 1 it connects with varying PEc/p classical 
thermal activation with MQT [fig. El arrow (b)]. Apart from limiting cases, pn|) is easily 
evaluated numerically; some results are shown in fig. |3] Hence, the QSE approach allows to 
give, in a seemingly transparent manner, a complete description throughout a broad range 
in parameter space and must be supplemented only for very low temperatures (T 0) by 
more sophisticated techniques jscJ ). 



C. Electron transfer in condensed phase 



The most prominent example of molecular electron transfer (ET) is the primary step 
of photosynthesis, but it is of fundamental relevance in a large variety of chemical and 



biological compounds 



35 1 . In the simplest case, two electronic sites denoted as donor and 

n 

acceptor, respectively, interact with a bath of harmonic oscillators (spin-boson model |36|). 
Completely equivalent is a formulation where the electronic system is coupled to a collective 
vibronic degree of freedom (Reaction Coordinate, RC) embedded in a thermal environment 
37[ , see fig. El The theoretical foundations to describe these systems have already been laid 
in the 1980s with the pioneering work by Marcus et al. 3^. It turned out that qualitatively 
two dynamical domains must be distinguished: For vibronic dynamics fast compared to 
the bare ET, characteristic for long-distance ET, one speaks of nonadiabatic reactions; in 
the opposite case of very slow vibronic motion, realized e.g. in mixed-valence compounds, 
the reaction is said to be adiabatic. Accordingly, for a donor-acceptor system transfer 
rates have been successfully calculated in the nonadiabatic regime by invoking golden rule 
techniques and for adiabatic reactions at sufficiently high temperatures by Kramers' flux 
over population method These quite different approaches reflect the different physical 
processes that control the transfer, namely, in the former case the electronic coupling between 
diabatic donor /acceptor states, in the latter one the sluggish activated bath motion on the 
lower adiabatic surface. For sufficiently high temperatures and slow to moderately slow 
bath modes, first Zusman 0] and later Garg |^ derived equations of motion for the 
electronic dynamics interacting with a damped RC treated in the classical Smoluchowski 
limit. Since these Zusman Equations (ZE) allow to derive a rate expression which to some 
extent comprises adiabatic as well as nonadiabatic effects, it is important to look for a low 
temperature extension Q]. 

The corresponding Hamiltonian is of the form with a system part Hs = H^i^ + i^Rc 
where 

Hel = -y^^ - -Y^'^ 

represents the bare electronic two-state system (EL) as an artificial spin-| system with donor 
|— ) and acceptor |+), a bias he, and an electron transfer coupling ^A; the second part 

Hrc = ^ + - coga. (45) 

describes a harmonic collective vibronic degree of freedom linearly coupled to the electronic 
system and to a heat bath (see fig. What one is interested in is the dynamics of a state 








q 



FIG. 5: Diabatic potential surfaces for the reaction coordinate (RC) in a ET process. A wave 
packet is excited by a laser pulse from a dark state (thin line) to the donor state (thick line, 
minimum at —qo). There it evolves on the coupled donor and acceptor (thick line, minimum at 
+qo) surfaces where the ET occurs in the Landau-Zener region near q = 0. 

initially prepared on the donor surface. It is given by the exact path integral expression for 
the reduced density matrix of the EL+RC system pafsil, q' a,f3 E { — ,+} |4Q|. 

Now, for strong damping (slow RC) a simplification according to the scenario discussed 
above is applicable. The difficulty then lies in the coupling between overdamped RC and 
electronic motion which basically makes a transition from one electronic state to the other 
only possible in the Landau-Zener (LZ) region where the harmonic surfaces of the RC in- 
tersect (see fig. El)- Thus, three additional time scales enter the game: The typical time for 
thermal activation into the LZ range, tact, the typical time the RC diffuses, after activation, 
in the LZ range, /^lz, and the typical time the density matrix resides in a non-diagonal state, 
tbiip- It turns out that (fTT|l must be extended to 



in order to reduce the path integral expression. This way, one derives for the probabili- 
ties Pai3{q,t) = Paf3{(l,(l,t) the extension of the ZE to lower temperatures, the Generalized 
Zusman Equations (GZE), 




1 h/3 



(46) 



£^P^^(g,t)±^ [P_+(g,t) 
P^±(g,t)±^[P-(?,t) 



P+-iQ,t)] 




(47) 



Here, the quantum Smoluchowski operators [see ()27|)] read 

1 d 



Mujl{q-7]qo) + - — 



(48) 



M-fdq 

and describe the overdamped quantum dynamics of the harmonic RC on the donor (77 = — ), 
the acceptor {rj = +), and the averaged potential surface {rj = 0). Note that the effect 
of quantum fluctuations also shows up in the Cq/k dependent coupling terms of the non- 
diagonal matrix elements. The coefficient n contains the equilibrium variance of a damped 
harmonic oscillator 

t^=-l^^l + Mulf3A. (49) 

By invoking again the tinae scale separation (j46p the electron transfer rate can be de- 
rived from the GZE 0, fl 
ward-|-backward) rate 

k 



42]. In case of a vanishing bias one finds for the total (for- 
A2 



l+g 

with an adiabaticity parameter 



exp {-(3E,/Ak) (50) 



g = TTK ^— — 51 

and with the reorganization energy E,. = 2cq/MuJq. For g ^ 1 one recovers the adiabatic, 
for (7 <^ 1 the nonadiabatic rate constant. The above rate expression looks like the classical 
Marcus/ Zusman result with a renormalized reorganization energy E^ —>■ Ej./n. Note that this 
simple renormalization only appears if quantum fluctuations in the Co-coupling terms of the 
GZE are properly taken into account. Since n > 1 quantum fluctuations always reduce the 
effective energy barrier which is to be surmounted. Of course, in the high temperature limit 
(k 1) we regain the known result, for lower temperature, however, significant deviations 
are observed, (i) The ratio k grows with decreasing temperature, at very low temperatures 
roughly linearly with f3. As a consequence, keeping all other parameters fixed, g becomes 
larger with lower T meaning that one approaches the range where the transfer is dominated 
by adiabatic processes at smaller values of A compared to the classical range, (ii) At lower 
T, the exponent in the rate expression f^Ej-/ n tends to become temperature independent in 
contrast to the classical result. Since k enters the exponent in ()5U|) . even relatively small 
deviations from the classical behavior K,/f3 = 1/(3 substantially influence the rate. In fact, 
if formally the limit of very low temperatures 7^/? ^ l,ujQh(3 ^ 7/1^0 is taken [which in a 
strict sense is out of the range of validity of (I^Uj)]. k/ [3 saturates and one obtains 

(3E, TT-f E, 



4ft 8u;o In(7/cuo) Tiuiq 



(52) 
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FIG. 6: Electron transfer rates vs. temperature for a symmetric system (e = 0) according to the 
expression (|50j) (solid line), the classical result {k = 1, dashed line), and precise Quantum Monte 
Carlo data (diamonds, from 0). Er/nA = 10 and A/{uj^/'j) = 2. 

which is identical to the overlap of two Gaussian wave packets with overdamped harmonic 
variance localized around ^Cq/MuJq. Hence, nuclear tunneling is included in the above rate 
formula, at least in a nondynamical way. In comparison with precise path integral Monte 

n 

Carlo results (from [43]) a remarkable agreement is seen over a broad temperature range, 
fig. El At low T, the rate enhancement due to nuclear tunneling is substantial. 



IV. CONCLUSIONS 



The description of quantum Brownian motion is still a challenging problem. In general 
it does not allow for "simple" solutions in terms of tractable time evolution equations and 
even a direct numerical evaluation of the exact path integral expression for the reduced 
dynamics is not always feasible. While in the limit of weak friction various master equa- 
tions have successfully been applied in the past, the opposite limit of strong friction has 
been left basically untouched. However, the naive expectation is that a strongly condensed 
phase environment quenches quantum fiuctuations, thus forcing the system to behave more 
classically. In fact, this is only one facet of the problem since quantum mechanically small 
fiuctuations in position are immediately accompanied by large fiuctuations in the conjugate 
momentum. A systematic simplification within the path integral formulation leads indeed 
to an equation of motion which to leading order coincides with the classical Smoluchowski 
equation, but contains substantial quantum corrections. 



Applications for chemical, mesoscopic, and soft matter systems have shown the decisive 
impact of corresponding quantum mechanical fluctuations on transport properties even up 
to high temperatures where thermal energies exceed relevant time scales of the damped 
dynamics. So far, this analysis has focused on one-dimensional systems, but by exploiting 
the developed techniques extensions to two and three dimensions are certainly feasible. 
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